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ABSTRACT 

General models describing the interactions between a pair of piezoceramic patches and 
elastic substructures consisting of a cylindrical shell, plate and beam are presented. In each 
case, the manner in which the patch loads enter both the strong and weak forms of the time- 
dependent structural equations of motion is described. Through force and moment balancing, 
these loads are then determined in terms of material properties of the patch and substructure 
(thickness, elastic properties, Poisson ratios), the geometry of the patch placement, and the 
voltages into the patches. In the case of the shell, the coupling between bending and in- 
plane deformations, which is due to the curvature, is retained. These models are sufficiently 
general to allow for potentially different patch voltages which implies that they can be suit- 
ably employed when using piezoceramic. patches for controlling system dynamics when both 
extensional and bending vibrations are present. 
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1 Introduction 


The use of piezoceramic elements as sensors and actuators has burgeoned in the last several 
years in applications ranging from the measurement and damping of vibrations in large flexi- 
ble structures to the control of noise in structural acoustics settings. Their utility as sensors 
derives from the property that when the element is subjected to a mechanical strain, a voltage 
proportional to the strain is produced. Conversely, they also exhibit the phenomenon that 
an applied polarization voltage across the unconstrained element produces in-plane mechan- 
ical strains in the material. Because of these properties, piezoceramic elements have found 
increasing success both as sensors such as strain gauges and accelerometers and as distributed 
actuators. Their success as actuators is augmented by the fact that they can be used to di- 
rectly control local vibrations without applying rigid body forces and torques, and due to their 
distributed nature, they are less prone to spillover effects in many control strategies. More- 
over, the piezoceramic elements or patches are inexpensive, lightweight, space efficient and 
can be easily shaped or bonded to a variety of surfaces. Hence a large number of patches can 
be used to sense and control without significantly changing the mass or dynamic properties 
of the system. 

In order to obtain optimal results with the piezoceramic elements or patches in sensing 
and control applications, it is necessary to have accurate models of the mechanics of induced 
strain actuation. This modeling also provides knowledge of the physical limitations of the 
piezoceramic patches as actuators in various settings. Detailed models have been developed 
for piezoceramic patch interaction with Euler-Bernoulli beams [2, 3, 4, 5, 10] and thin plates 
[6, 12]. Because many of the initial applications of piezoceramic elements were in settings 
involving the sensing and control of bending deformations (these vibrations are dominant in 
many low frequency vibration and noise control problems), most of these models concentrate 
on patch configurations which excite pure bending motion of the substructure with more 
limited discussions of pure extensional excitation. It was not until [10] that a model was 
developed which provided for simultaneous excitation of both bending and extensional de- 
formations in an Euler-Bernoulli beam. One motivation for developing such a model is the 
observation that in complex coupled systems, in-plane vibrations with small displacements can 
have large in-plane energy levels due to the property that beams are much stiffer in extension 
than in bending. This in-plane energy can then couple into flexural vibrations at structural 
discontinuities such as joints, thus necessitating the control of both bending and extensional 
vibrations in such structures. As determined by Fuller et al [8] through experimental work, 
simultaneous reductions in both flexural and extensional deformations in a beam can be ob- 
tained through the use of asymmetric pairs of piezoceramic actuators and sensors in adaptive 
control schemes, and the analytic work in [10] was a first step toward developing a model 
which could be used in further such control settings. In that work, force and moment balanc- 
ing were used to determine expressions for the moments and strains induced by the activation 
of a single piezoceramic patch which was bonded to an Euler-Bernoulli beam. 


In addition to beams and plates, thin elastic shells are often used to describe various 
structural components as well as when modeling the coupling between structural vibrations 
and their radiating or receiving acoustic fields. For example, the transmission of sound through 
an airplane fuselage due to low frequency, high amplitude exterior acoustic fields can be 
modeled by a vibrating thin cylindrical shell which is coupled to an interior acoustic pressure 
field [9]. In order to optimally control the interior noise via piezoceramic patch actuation, 
one first needs to accurately model the interactions between the patches and the shell. This 
raises modeling issues which differ from those encountered in the beam and plate analyses in 
that the in-plane and bending vibrations are coupled in the cylindrical shell due to curvature 
effects. 

Analytical models describing piezoceramic. patch/cylindrical shell interactions have primar- 
ily been based on layered shell theory [11, 19] or the use of flat plate piezoceramic coupling 
results when determining the resulting loading on the shell [15]. In the first case it is assumed 
that the piezoceramic material makes up an entire layer of the elastic structure and hence 
this model is of limited use when considering small patches as actuators. When using the 
flat plate theory, it is assumed that the patch dimensions are small in comparison with the 
cylinder radius. Curvature properties are then neglected when modeling the coupling between 
the patch and shell and determining the loading due to activation of the patch. 

In this work, we present general models for the interactions between a pair of piezoceramic 
patches and elastic substructures consisting of a beam, plate or thin shell. In the case of a 
shell, the patches are assumed to be curved and the coupling between bending and in-plane 
deformations, which is due to the curvature, is retained. The techniques used to develop the 
shell/patch interaction model are also used to develop general models describing the moments 
and forces which are generated by the activation of piezoceramic patches which have been 
bonded to a flat plate or beam. These models differ from those in [6, 12] and [2, 3, 4, 5] in 
that they allow for different voltages into the individual patches thus admitting the analysis 
of simultaneous excitations of both bending and extensional components. In the case of the 
beam, the model is slightly more general than that in [10] since it is derived for two active 
patches. Hence it is appropriate for models in which both patches actuate with potentially 
different voltages into the patches. In the case of one actuating patch, how 7 ever, the model 
reduces to that in [10]. From a control perspective, these models are important since they 
provide for greater latitude in designing control strategies involving the use of piezoceramic 
elements to affect both the bending and extensional properties of a structure. 

As a prelude to the development of the patch interaction models, equations of motion for 
the underlying substructures are presented with special attention paid to the contributions 
due to externally applied moments and forces since this is where the interactions between 
the patches and substructure occur. The analysis leading up to the structural equations also 
motivates many of the techniques which are used to develop the patch interaction models. 

To this end, a synopsis of the derivation of the strong form of the time-dependent Donnell- 
Mushtari thin shell equations from Newtonian principles (force and moment balancing) is 
presented in Section 2. A complete treatment of this topic can be found in [13, 14, 16, 17, 
18] and the included discussion is limited to summarizing that material which is needed for 
developing the shell/patch interaction model as presented in the following section. The choice 
of the Donnell-Mushtari model is for ease of presentation and as noted at various points in the 
discussion, the patch /shell interaction model can be easily extended to higher order models 
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as warranted by the physical situation. 

An inherent disadvantage of the strong form for the equilibrium equations when the ex- 
ternal loads are generated by piezoceramic elements is the resulting presence of the first and 
second derivatives of the Heaviside function due to the finite support of the patches. As a 
result of this as well as other identification and approximation issues, we then develop the 
weak form of the time-dependent Donnell-Mushtari shell equations. This is done in more 
detail since this development is less readily available in the the literature. This formulation is 
advantageous in many approximation schemes, admits the identification of discontinuous ma- 
terial parameters, and eliminates the problem of differentiating the Heaviside function since 
the derivatives are transferred onto the test functions. 

The second section concludes with a synopsis of the strong and weak forms of the Kirchhoff 
plate and Euler-Bernoulli beam equations. As in the shell discussion, particular emphasis is 
placed on the contributions of externally applied forces and moments since this is where the 
coupling between the substructure and piezoceramic patches occurs. 

The patch contributions to the cylindrical shell equations are developed in Section 3. In 
order to determine the loading due to patch moments and forces, it is useful to first express 
them in terms of the normal strains and changes in curvature of the middle surface of the 
cylindrical shell. To do this, the stress-strain relations in the patch and shell are developed 
which then allows the moment and force resultants for the patch to be formulated in terms of 
midsurface shell properties. The unknown normal strains and midsurface changes in curvature, 
and hence the patch moments and forces, are then determined by moment and force balancing. 
In this manner, the loading due to activation of the patches can be expressed in terms of 
material properties of the patches and shell (thickness, elastic properties, and Poisson ratios), 
the radius of curvature of the shell, and the voltage being applied to the patches. 

In Section 4, the techniques of the third section are tailored to composite structures consist- 
ing of piezoceramic patches which are bonded to plates and beams. The resulting plate/patch 
interaction model is shown to be equivalent to that of [12] in the special case when pure 
bending motion is excited (the model in [12] was derived by isolating the interface stress for 
the system and treating it as the unknown to be determined). Due to its generality however, 
our model also allows for more complex interactions involving both bending and extensional 
components since the voltages into the individual patches can differ. As discussed earlier, the 
beam/patch interaction model reduces to the model in [10] in the case of one actuating patch 
but is slightly more general in that it also admits models in which two patches are used for 
actuating with potentially different voltages into the patches. As with the shells, this provides 
structure/patch interaction models which can be used in various structural and structural 
acoustics control settings. 
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2 Underlying Shell, Plate and Beam Equations 


Throughout this discussion, we consider a thin circular cylindrical shell of radius /?, thickness h 
and having the axial coordinate x as shown in Figure 1. The variable 2 measures the distance 
of a point on the shell from the corresponding point on the middle surface (z = 0) along the 
normal to the middle surface. 



Figure 1 . The Cylindrical Thin Shell. 


Strain-Displacement Relations 

By combining Love’s shell assumptions with the strain-displacement equations of three 
dimensional elasticity theory, one obtains the strain relations 


C# — T % K j 

1 


ee = 


7x0 


1 + z/R 


(e 0 + zk b ) 


1 +z/R L 


£ x e + 2(1 + 


( 2 . 1 ) 


where e x and e 0 are normal strains at an arbitrary point within the cylindrical shell and 7^ is 
the shear strain. Here e x , e$ and e x0 are the normal and shear strains in the middle surface and 
k x , k 0 and r are the midsurface changes in curvature and midsurface twist (see [ 14 ], page 8). 

Note that within the framework of infinitesimal elasticity, the equations (2.1) are exact 
and in the Byrne-Fh’igge-Lur’ye shell theory, these represent the exact form of the kinematic 
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equations. In the Donnell-Mushtari theory, one neglects the underlined terms zj R with respect 
to unity thus leaving 

C x £ x "T ZK X 

ze — + zkq ( 2 . 2 ) 

7xs = e xe + z (l + ^) r . 

In terms of the axial, tangential and radial displacements u, v and in, respectively, the 
expressions for the midsurface strains and changes in curvature for the cylindrical shell are 


du 



d 2 w 

dx 


dx 2 

1 dv w 

UTe + r 

, K-0 

1 d 2 w 

R 2 d6 2 

dv 
d x 

1 du 


2 d 2 w 

R~d0 

> T = ” 

RdxdO 


1 dv 

+ jrM 

2 dv 
R dx 


(2.3) 


As before, the underlined terms are retained in the Byrne, Fliigge and Lur’ye theory 
and are discarded in the Donnell-Mushtari theory. We point out that the equations (2.1) 
and (2.3) differ from those arising in the theory of flat plates both in the presence of the 
length differential RdO as well as in the retention of the strain terms e x and e# (only bending 
contributions are considered in the corresponding models of the transverse vibrations of a flat 
plate). 


Stress-Strain Relations 


To determine the constitutive properties of the shell, it is assumed that the shell material is 
elastic and isotropic. Hooke’s law in conjunction with the assumption that the transverse shear 
stresses a xz and cr$ z as well as the normal strain component e z are small in comparison with 
other stresses and strains (these conditions are part of Love’s third and fourth assumptions) 
then yields 

E , 

cr x = 7 y( e * + veo) 

1 — v L 

E 

ere = z(e e + ve x ) (2.4) 

1 — v l 

ax& = <T0 * = 2 (l + «/)' 1x6 

where a x and a$ are normal stresses and <j x q and a$ x are tangential shear stresses. The 
constants E and v are the Young’s modulus and Poisson ratio for the shell. 
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Force and Moment Resultants 


By integrating the stresses over the face of a fundamental element, the force resultants 
be expressed as 


' N x ' 

N xQ 

. Q* 



<7 x 
<7x0 
<7 x* z 



dz 


can 


and 

’ No 
Nox 

_ Qe 

Similarly, the moment resultants are 



<70 

<7 Ox 

<7Qz 


dz . 


(‘ + i) zdz 


and 

' Me 
Mg x 

The orientation of the various forces and moments are shown in Figure 2. We point out that 
the transverse shear stresses a x ~ and ag z are used when obtaining the force resultants Q x and 
Qe even though they are omitted in the constitutive relations. This is one of the contradictions 
which arises in the classical shell theory. 


rh/2 

<70 

J —h/2 

<7 Ox 


M x 

= /‘ /2 ^ 

Mxo _ 

J—h/2 <7x0 


Qx 



Figure 2. Force and Moment Resultants for the Cylindrical Shell. 
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In the Donnell-Mushtari theory, the underlined terms z/R are neglected in comparison to 
unity, and the integrals are determined accordingly to yield 


T Eh du ( 1 dv w'' 

Nx = (1 -v 2 ) \dx + " \R99 + ~Rj 


N 0 = 


Eh 


(1 - /+) 

N x e — Nq x = 


1 dv w du 

1 b V 

R 86 R dx 


Eh 


2(1 + v) 


dv 1 du 
d^ + ~R99 


M x = 


M e = 


Eh 3 


9 2 w v d 2 w 


+ 


12(1 -v 2 ) [dx 2 R 2 d6 2 

Eh 3 1 8 2 w d 2 w 
12(1 -v 2 ) [R 2 'W + V ~d^ 


M x g = M 0X 


Eh 3 


d 2 


w 


12/2(1 + v) 8x96 


(2.5) 


Similar expressions are obtained in the higher order theories. 


Strong Form of the Donnell-Mushtari Shell Equations 


The equations of the dynamic equilibrium of the element are obtained by balancing the 
internal force and moment resultants as shown in Figure 2 with any externally applied forces 
and moments. Let 

<7 — Qx^'x T (]9^9 T (jiihi 

and 

m = m x i x + mete 

denote the surface forces and moments due to an external field which is acting on the middle 
surface. Hence q and m have units of force and moment per unit area, respectively. 
Considering equilibrium of the forces in the x , 9 and z directions yields 


n dN x , dN ex , 

R~ — I — ™ h R<Ix 


dx 


80 


0, 


dN e + R N* + Q $ + R j 0=:Oi 


R 


99 

9Qa 


dx 


+ 


dx 

dQe 

99 


No + R(jn — 0 , 


( 2 . 6 ) 


respectively. In the Donnell-Mushtari theory, the transverse shearing force Qq is considered to 
be negligible in the second equation of (2.6) and is subsequently neglected when determining 
the final equilibrium equations. Similarly, with 0 as a reference origin, the balancing of 
moments with respect to 0, x and z yields 


D dM x _ 9M gx , D ; _ n 

R T r\/\ RQ x T Rm$ — 0 , 


dx 


99 


9 M e D 9M x6 , D . _ n 

— b R — x RQe + Rm x — 0 , 

o9 ox 

N x g - Ng x - ^ = 0 , 


(2.7) 
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respectively. By referring to the integral definitions of N x $, N$ x and Mg x , it can be seen that 
the third expression in (2.7) is identically satisfied due to the symmetry of the stress tensor. 

Time enters the equilibrium equations through the inertial terms; hence for time-dependent 
problems, the force q x is replaced by — ph%^-\-q x where p is the density in mass per unit volume 
of the shell. Similar substitutions are made for q$ and q n . By combining (2.6) and (2.7), one 
arrives at the time-dependent Donnell-Mushtari equilibrium equations for a thin cylindrical 
shell with radius of curvature R and thickness h 


R ph ^-R™- 


dN, 


0x 


Rph 


dt 2 

d 2 v 


dx 

dN e 


dO 


Rtfx 


dt 2 dO 


R^- = Rqo 

OX 


(2,8) 


Rph 


d 2 w 

dt 2 


- R 


0 2 M X 

dx 2 


1 d 2 M e o d 2 M x8 

r do 2 1 dxdo 


T Ng — Rq n T R 


drhg 

dx 


+ 


dm x 

~W ' 


We note that the representation of the external loads as surface moments and forces is 
convenient when deriving the strong form of the equations of motion. However, in many 
applications where it is necessary to actually determine expressions for these loads or when 
using the weak form of the equations, it is advantageous to represent these loads in terms 
of line forces and moments. To accomplish this, let M x , Mg, N x , and Ng denote the external 
resultants acting on the edge of an infinitesimal element which have the same orientation as 
the internal resultants depicted in Figure 2 (with units of moment and force per unit length 
of middle surface). Force and moment balancing can be used to write the area moments and 
in-plane forces in terms of these line moments and forces, thus yielding 


(jx = 


rn x = 


di W x 

dx 

1 dMg 
R~df 


qe - 


mg 


1 dN e 
R~dO 
dM x 
dx 


(2.9) 


We point out that the first expression in (2.9) can be obtained from (2.6) simply by 
replacing N x by N x and deleting in the first expression of (2.6). Similar analysis leads to 
the other expressions in (2.9). The use of these line moments and forces in (2.8) is equivalent 
to including the external resultants directly when determining the equations of moment and 
force equilibrium for an infinitesimal shell element as done in (2.6) and (2.7). 

The substitution of the internal moments and forces in (2.5) and the external resultants 
from (2.9) then yields 


1 d 2 u d 2 u 1 — u d 2 u 1 + v d 2 v v dw (1 — u 2 )dN x 
Cl~di 2 ~d^~ 2 R 2 W 2 ~ 2 R dxdO ~ ~RlN = Eh d7 


d 2 i 


1 


d\ 


1 d 2 v 1 + v d 2 u 


a a? 


2 i)x'' R?9P 2 R dxdO 


1 dw 

lR~dO 


d 2 w v du 1 dv 1 h 2 4 (1 — v 2 ) 

(- 1 b — w A Vue = 

L uo- Rdx R 2 dO R 2 12 Eh 


cl dt 2 


qn 


( 1 — v 2 ) 1 dN e 
Eh R~df 

1 d 2 Mg dM x 


( 2 . 10 ) 


R 2 do 2 


dx 2 
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where again, u, v and w are the axial, tangential and radial displacements, respectively [14]. 
The constant Cl given by 

r = h 
L Ua -"?) 

is the phase speed of axial waves in the cylinder wall. The external line forces N x and N$ 

A A 

and moments M x and M$ have units of force and moment per unit length of middle surface, 
respectively, and are generated in our problem by the activation of the piezoceramic patches. 
The load q n is left as a surface force since this is the form that it usually takes in problems 
involving the excitation of a shell through normal forces (an example of a normal force in this 
form is the pressure exerted on the shell due to an exterior or interior acoustic field). 

We again emphasize that the resultant expressions in (2.9) (and hence the system (2.10)) 
were derived for an infinitesimal element; hence certain modifications must be made when 
considering the global form of the resultants and equations (as is necessary when the resultants 
are generated by a piezoceramic patch). In certain cases (e.g., for certain types of moments 
and forces), the system (2.10) agrees with the strong form of the global shell equations. In 
general, however, this is not true, and one must exercise extreme care in determining the form 
of the global representations for the moments and forces. 



Weak Form of the Donnell-Mushtari Cylindrical Shell Equations 


In order to find the weak form of the shell equations, the kinetic and strain energies of 
the shell are needed. By combining the Kirchhoff shell hypothesis with the strain results 
from classical elasticity theory, it follows that the strain energy stored in the shell during 
deformation is given by 


\ thj 2 r2ir rl 

o/ / / {Cxe x + ageg + a x 0-f x e) {\ +z/R) 

Z J-h/2 j 0 Jo 


RdxdOdz 


where the stains and stresses are given in (2.1) and (2.4), respectively. Substitution and 
integration (with (1 +^/i?) _1 replaced by its geometric series expansion and neglecting powers 
of z in the integrand which are greater than two) yields 




h? \ 2 ( r 1 2 \ 2 

T j ^ (^x "h 2(1 Z'') I k x Kq ” 1 ~b x^x q ) 

_<i zAhi r + 4 + \fldxde . 

2 R R 2 2 R 2 I 


With the change of variables s = x/R ) the total strain energy can be written as 


1 r 27r R! R Eh 

2 Jo Jo (1 — is 2 ) 


[Idm + hi bfl\ dsdd 
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where k = h 2 /(V 2 R 2 ), lp\i is the integrand corresponding to the Donnell-Mushtari theory 
and Ibfl denotes the terms which are retained to yield the Byrne, Fliigge and Lur’ye strain 
energy. These two components are given by 


, fdu dv V on \ 


du 

ds 


+ k | ^V 2 te) 2 — 2(1 — u) 


d 2 w d 2 iv 
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i / dv < 9 u \ 2 
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du d 2 w 


1 H d 9 j ^ ds ds 2 


+ 2 w 


d 2 


w 


de 2 


+ U) 


For simplicity of presentation, a weak form of the shell equations will be developed using 
the Donnell-Mushtari strain expression; a corresponding set of equations can be derived in a 
similar manner in the Byrne, Fliigge and Lur’ye case. 

The kinetic energy of the shell is given by 





R 2 dsd 9 . 


Throughout this development, it is assumed that the shell satisfies shear diaphragm bound- 
ary conditions at x = 0 ,/; that is, it is assumed that 


v — w = N x — M x = 0 


at the ends. This is done merely to demonstrate the equivalence between the weak form which 
follows and the strong form already discussed; other boundary conditions can be treated with 
similar arguments. It should be noted that the conditions v = w = 0 at the ends are essential 
boundary conditions and hence must be enforced on the chosen state space. 

For an arbitrary time interval [/o,/i], consider the action integral 

A\u] = [ ' (T - U)dt 
Jto 

where u — [u } v y w] is considered in the space V = fl) x Hl(Q) x H 2 {yt). Here fi denotes 
the shell and the subscript b denotes the set of functions satisfying the essential boundary 
conditions. One then considers variations of the form 


( 2 . 11 ) 


u(t , 7 ', 0 , x) 


’ Vi ' 

v(t, 7 ’, 0 , x) 

+ £ 

772(0^2(7', 0 ,x) 

w(t , 7 ’, 0 , X ) 


_ 773(0^3(7-, 0 , x) _ 


Here ff = [771,772,773] and <f) = [<^i, fa, ^3] are chosen so that 
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i. ) u(t, •, •, •) € V 

ii. ) u(to, ■, •) = u{t\ • 

Note that this enforces if £ [// 2 (0, 7 1 )] 3 , rffto) = if(t\) and <f £ V 7 . 

Hamilton’s principle states that the motion of the shell must give a stationary value to the 
action integral when compared to variations in the motion, thus leading to the requirement 
that for all <j>, 

-A [u + £$ 


de 


0 . 


£=0 


With the definition (2.11) for the action integral, Hamilton’s principle leads to the condition 




r2n ff/R 

= / Ph 

:=0 Jto JO JO 


du dift dv di ]2 , dw cb/ 3 

’9\ + 3 


dt dt 


dt dt 


dt dt 


u f 2n f (/R Eh 


ft 1 fZ 7T f 

J \q Jo Jo 


( du dv 


(i - u 2 ) \{ds + de + ( ?/1 ds + 7/2 de 


dd>. 


R 2 dsdO dt 

d<j> 2 


+ mfo 


-(l-n) 


d(f > t dv d<j > i du d<j )2 du 

mw 17 + m~§7 + r, W + m T,* 3 


l f dv du\ ( d<f >2 d<f> i\ 


+k< J7 3 V 2 wV 2 <^3 - (1 - v) 


d 2 w d 2 <f > 3 d 2 w d 2 4>3 

m W : ~d^ r + 


“ 27/3 


d 2 w d 2 (j >3 


dsdO dsdO 


>dsdOdt 


Note that this must hold for all arbitrary intervals [f 0 , ii] and all admissible perturbations. 
Temporal integration by parts in the first integral in conjunction with the underlying condition 
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that ij(t 0 ) = T]{ti) then yields the coupled system of equations 
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The weak form of the equations of motion for the unforced shell is thus 
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law* + U + "® + "“j + + - ,) (* + « j ¥|" 0 


f 27r fdfi f fi 2 d 2 u / clu d» \ d <^2 U. _ H — , 

l>2 + + 7T7 + ^ a/l + v ) + 


I f 


dv du \ d(f>2 


Cj dt 2 


ds dO 


de 


ds de ds 


dsdO = 0 


2* f t/R f Q2 W ( du dv \ . 

<h+ a2 + w ) & 


If 


a dt 2 


ds de 


+k^ 2 wV 2 cf> 3 -(1 - v) 


d 2 wd 2 <f> 3 , d 2 wd 2 <j > 3 0 d 2 iv d 2 <t>3 

d - ~ — ~ „ 


de 2 ds 2 ds 2 dO 2 


‘ dsde dsde 


}dsde = 0 


I j 2 

for all <j> G V . Again, the constant Cl = is the phase speed of axial waves in the 

cylinder wall. 

In terms of the moment and force resultants (see ( 2 . 5 )) and the original axial variable x, 
the weak form is 

r2v ft ( d 2 U , . 3<f)\ 




f [ { Rpk §£* + N, d -t + RN„^ i dxdO = 0 


dxdO = 0 


d<j) 2 


de 


f 2 v ft ( d 2 w 

L L Ri,h w h + N,h ~ 


dx 

1 Uxde = °- 


(2.12) 


d 2 ^»3 

’dxdO 
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The derivation thus far has been for the unforced shell. To include the contributions of 
applied external forces and moments which do nonconservative work on the shell, one can ap- 
peal to an extended form of Hamilton’s principle or more formally include these contributions 
directly in the system (2.12). Both techniques yield identical final equations and for ease of 
presentation, we will take the latter approach. 

The inclusion of the applied line forces and moments N x , N$ an d M x , M$ and the surface 
load q n in the system then yields 




d \2 + N e ^- + RN x9 


86 

d<t> 2 


RN x ^-\ dxde = 0 


n> 


8 2 


w 


8t 2 


86 


<f> 3 + Neh — RM. 


dx 

8 2 fa 


- Nt 


x dx 

d(f> 2 


} 


1 


86 

8 2 fa 


dxdO = 0 


* dx 2 R Me 86 2 


— 2M x g 


ddfo 

8x86 


-RM 3 + + ^M,^\dxde = o 


(2.13) 


for all <£ € V as the weak form of the Donnell-Mushtari equations of motion for the forced 
she 11 - 

With the assumption of sufficient smoothness, the weak solution in this form is consistent 
with the strong solution in (2.8). The vanishing of several of the boundary terms which arise 
during integration by parts is a result of the choice V = II l (ff) x Hi (fl ) x H 2 ($l) for the 
function space since the state variables and test functions are required to satisfy the essential 
boundary conditions 

v — w = 0 

at x = 0, 1. 

We point out that in the weak form (2.13), one is not required to differentiate the applied 
force and moment resultants N x , N$, M x and Mq as is required in the strong form (2.10). This 
proves to be very beneficial when these terms are generated by finite piezoceramic patches as 
discussed in the next section. 


Plate Equations 

Consider a thin rectangular plate whose edges lie along the coordinate lines x = 0 , 1 and 
y = 0, a. We assume that the plate subjected to both longitudinal and transverse loading 
via the surface forces and moments q x ,qe,q n and rh x ,rhg. With u, v and w denoting the 
displacements in the x,y and normal directions, respectively, the strong form of the Kirchhoff 
plate equations is given by 
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ph 

ph 


d 2 u dN x ON , 


yx 


dt 2 

d 2 v 

dt 2 


dx 


dy 


— <7x 


dN v N s 


xy 


d 2 w 

P k ~^2 ~ 


dy 

d 2 M x 


dx 


~~ Q y 


(2.14) 


d 2 m„ d 2 M. 


xy 


d 2 M 


yx 


dt 2 dx 2 dy 2 dxdy 

where the moment and force resultants are 

Eh 


dydx 


— (jn + 


dm x dm ,, 


dy 


dx 


Nr 


— v* 


Eh 


M — 

V 1 - 1/3 


f du dxf 
K dx + 1 dy , 

'dv dxE 

Aj + 


N 3 


xy 


N, 


Eh 


yx 


2(1 +v) 


dv dvr 
dx + dy. 


M x = ■ 

Eh 3 / d 2 w 

12(1 — u 2 ) \dx 2 

My = - 

Eh 3 (d 2 w 

12(1 - v 2 ) \ dy 2 

II 

„ Eh 3 

Myx ~ 12(1+// 


+ v 


d 2 


w 


+ V 


dy 2 , 

d 2 w' 


dx 2 

d 2 w 


The first two equations in (2.14) describe the longitudinal movement of the plate while the 
third equation describes the transverse motion of the plate. 

To find the weak form of the equations, the vector u — [u, iv] containing the displace- 
ments in the x, y and normal directions is considered in the space V = Hl(0) x HUa)xHm 
where f l denotes the plate and the subscript b denotes the set of functions satisfying essential 
boundary conditions for a specific problem. By using analysis similar to that just described 
for cylindrical shells, the weak form of the equations of motion for the plate can be found to 

' yx dy 


a: 

ff 

n‘ 


f „ 3*1 , Oh fy a* 

ph m , *1 + N, g x + -V g a N x 8x 


dxdy — 0 


,S 2 V dh 34 

Ph — fo + 


2 _ * 3h 


’ dy 


’ l ^A. 

P h W h 


M ^2 
1 dx 2 




d 2 <f> 3 

dy 2 


j dxdy = 0 

d 2 (j> 3 


(2.15) 


* d 2 (b 3 - d 2 <b 3 

— 3 + M x Q + M y 


dx 2 


dy 2 


M i r ~ J _ m 

xy dxdy yx dxdy 


>dxdy — 0 


for all (j) = [<f>u<f>2i<f>3\ ^ V\ As in the case of the thin shell, the external line forces and 
moments N x , N y} M x and M y are related in an infinitesimal sense to the corresponding area 
forces and moments q x ,q y ,m y and m x appearing in the strong form of the equations by the 
relations 


Qx 


d_ N* 

dx 


% - 


m r = — - 


dMy 

dy 


.dNy 

dy 

dMr 

dx 


(2.16) 
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If the solution has sufficient smoothness, integration by parts can be used to show that the 
weak solution is consistent with the strong solution in (2.14). 


Beam Equations 


The motion of an undamped thin beam of length £ and width 1 can be determined from 
the dynamics of thin plate theory by considering only the vibrations in the x-direction along 
with the usual transverse vibrations (in the z direction). From (2.14) this yields the strong 
form of the Euler-Bernoulli beam equations 


where 


ph 

ph 


d 2 u 

dN x _ 

dt 2 

dx 

d 2 w 

d' 2 M x 

dt 2 

dx 2 

= Ekp 

OX 


Eh 3 d 2 w 


12 dx 2 


q n + 


dlhy 

dx 


= -El 


d 2 w 

dx 2 


(2.17) 


Note that / = h 3 / 12 is the moment of inertia for a beam of width 1. 

A corresponding weak or variational form of the equations can be determined by choosing 
V = Hi (0) x H 2 (Q) for the space of trial functions where 0 denotes the beam and the subscript 
b again denotes the set of functions which must satisfy the essential boundary conditions. 
Through either an energy derivation such as that given for the thin shell, or simply integration 
by parts, one arrives at the variational form 


jf + **17 “ dx = 0 for a " h s H ‘ (n) 

° 1 * 1 <2J8) 

£ - <j n fa + * = 0 for a" h € Hl(V) 

of the beam equations. We point out that in this form, one is not required to differentiate the 
external force or moment resultants, N x and M x , which proves to be very useful when these 
terms are generated by the activation of finite piezoceramic patches. 


3 Patch Contributions to the Shell Equations 


For a thin cylindrical shell, the strong and weak forms of the equations of motion are given 
by (2.8) and (2.13), respectively. In the case of the weak form, it is seen that the loads 
can be written in terms of the line forces and moments N x , N$, M x and Mg and the normal 
surface load q n , while the strong form contains surface loads and the derivatives of surface 
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moments. In the problem under consideration, these quantities result from the activation of 
piezoceramic patches of thickness T which are assumed to be perfectly bonded to a cylindrical 
shell of thickness h with midsurface radius R (see Figure 3). As shown in Figure 4, the patches 
are assumed to be situated so that their edges are parallel to lines of constant x and 9. Because 
the patches generate no shear strains, the exterior load q n is taken to be q n — 0. If the weak 
form (2.13) is used, the external line moments and forces are simply 


(3.1) 


M x — (Mx)pe , Me — {M#) pe 

N x = (N x ) pe , Ne = (A^)pe 

where ( M x ) pe , (Me) pe , (N x ) pe and (^) pe are the respective moments and in-plane forces which 
are generated by the patches. The subscript pe is used to denote patch properties and to 
help differentiate them from shell properties which have no subscript. When it is necessary 
to differentiate between the two patches, the outer will be denoted with a subscript pt\ with 
a subscript pt<i being used to denote the inner patch. 

However, if one is using the strong form (2.8) of the equations of motion with piezoceramic 
actuators, the surface moments and forces to be used in (2.8) are given by 

1 d(M e ) P e ^ _ d(M x ) pe 

(3.2) 


m T = 


R 86 


m 0 — 


dx 


q x — — Si } 2{ 1 * * * * * * * * X )S 1,2(6) 


d(N x ) 


pe 


qo = -Si, 2 (x)S h2 (6) 


1 d(N e ) pe 


dx ’ , ' 7 / p g e 

For a patch with uniform thickness and bounding values x x ,x 2 ,di and d 2 , the presence of the 
indicator function 

f 1 , x <(x x + x 2 )j2 

S\, 2 (x) = < 0 , x = (xi -f x 2 )j2 (3.3) 


1 -1 , x > (x\ + x 2 )/2 

derives from the fact that the forces generated by the patch in the x-direction are antisym- 

metric (equal in magnitude but opposite in signj about the line x = (jq + x 2 )/2. The same 

holds true for the forces in the 0-direction with §i i2 (6) being defined in an analogous manner 

to S , i i2 (*) in (3.3). 

We point out that the differences between the external surface force expressions in (2.9) 

and (3.2) are due to the fact that the former were derived for an infinitesimal element whereas 

the latter are global expressions which preserve the overall signs of the forces generated by the 
patches as well as reflect the discontinuities due to changes in sign. These differences result 
from the property that the sense of the forces is highly dependent on the specified location 
of the axis origin on the neutral surface. Hence the direction of forces throughout the patch 
differs in some locations from those observed in the infinitesimal element thus necessitating 
the inclusion of the indicator functions in (3.2). 

Unlike the forces, the action of the moments is specified with respect to a fixed point on 
the neutral surface (the point 0 for the element in Figure 2 or a point on the left edge of the 
shell in Figure 1). As long as the orientation of the infinitesimal element and full shell with 
patches are the same, the line moments derived for the infinitesimal element will be consistent 
with those of the full structure. Thus the expressions for the general infinitesimal moments 

in (2.9) need no modifications when describing the surface moments generated by the patches 
as given in (3.2). 
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Figure 3. Strain Distribution for the Composite Structure. 



In order to determine (M x ) pe , {Mg) pe , (N x ) pe and (Ng) pe and hence the loads on the shell, 
it is useful to write them in terms of the normal strains e x ,£g and midsurface changes in 
curvature k x , Kg of the middle surface (z = 0) of the cylindrical shell. That is, we want to 
express the patch moments and forces in terms of the reference surface characteristics of the 
cylindrical shell. 

We emphasize that due to the presence and activating nature of the patches when a voltage 
is applied, the normal strains and changes in curvature are no longer given by the expressions in 
(2.3) which were derived for a homogeneous thin cylindrical shell. At this point, e x ,£e, and 
Kg are considered to be unknown and are determined by formulating stress-strain relations in 
the patch and shell followed by the balancing of moments and forces in the combined structure. 
In this manner, expressions for these midsurface characteristics (and hence the resulting patch 
moments and forces) can be found in terms of the material properties of the shell and patch, 
the radius of curvature of the shell, and the voltages being applied to the patches. 
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Stress-Strain Relations in the Patch 


From (2.1), the infinitesimally exact strain relationships for a cylindrical shell with mid- 
surface radius R are given by 

(^x ~F ^^x) 

1 , , , (3-4) 

'* “ TT7JH (e " + ZK,} 

where e x and e$ are the normal strains at an arbitrary point on the cylindrical shell. If 
the patches and cylindrical shell are thin in comparison with the radius of curvature of the 
shell, then it is reasonable to assume that the relationship (3.4) is maintained throughout the 
combined thickness h + 2T as shown in Figure 3 (see also [12]). Hence we will take 


(^x)pe — (^-x + % K x ) 

, x 1 , , x (3.5) 

( e 0)pe = i+g J R ( £e + ZKe > 

where (e x ) pe and (c0) pe are the normal strains at an arbitrary point on the patch. Note that 
this assumption implies that the strains at the interface are continuous and that the centers 
for the radii of curvature for the shell and patch are concurrent. As seen in (3.4) and (3.5), 
the tangential strain distribution in the shell and patch is in general nonlinear in z (see also 
Figure 3). 

We point out that the model at this point differs from the flat plate piezoceramic coupling 
model [15] both in the presence of the term z/R in the tangential strain expression of (3.5) 
and in the fact that the model retains the coupling between the normal midsurface strains and 
the changes in curvature (this is analogous to simultaneously considering both longitudinal 
and transverse vibrations in a plate). Although the ratio z/R is neglected when deriving 
the Donnell-Mushtari model (see (2.2)), we retain it here so that curvature effects are fully 
included in the coupling between the patch and shell. The retention of this term also ensures 
that the patch interaction model can be directly applied to higher order shell models without 
necessitating changes to accommodate the greater accuracjc 

From the constitutive relations in (2.4), it can be seen that the stress distributions within 
the cylindrical shell are 

E 

= TTT 7 ( e * + ve o) 

E m 

oe = 7 o (eo + ve x ) . 

The stress distribution in the patches will contain contributions from both the free piezo- 
ceramic, actuator strain and the strain distribution in (3.5). At this point we assume that 
when voltage is applied and the patch is activated, in accordance with basic shell theory, 
equal strains are induced in the x and 9 directions and the radius of curvature is not changed 
in either direction. Patches satisfying this assumption could be made, for example, by taking 
a portion of a thin- walled tubular piezoceramic. element. For the outer patch, the magnitude 
of the induced free strains is then taken to be 


£pe i 


(f-x)pei — ( e e)pe i 
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where d^\ is a piezoceramic strain constant and V\ is the applied voltage. We point out 
that when a voltage is applied to a patch with edge coordinates X\,X 2 , 0 \ an( ^ ^ 2 ? the point 
(x,0) = ((xi + x 2 )/2,R(0i + 0 2 )/2) will not move whereas the axially symmetric points on 
either side will move an equal amount in opposite directions. This motivates the use of the 
indicator functions at various points throughout the development 

Assuming that the two patches have the same Young’s modulus, E pe , and Poisson’s ratio, 
z/ pe , the stress distribution in the outer patch is given by 

E 

(&x)pei “ “ ( e x T Vpe e 6 ~ (1 T Vpe) e pe ]) 

^ ^pe 

E 

{&$}pe i “ "j 2 (^0 T Vpe^x (^T Vpe^Zpe i) 

* ^pc 

with the negative signs resulting from conservation of forces. Similar expressions are used for 
the induced free strain and stress distribution in the inner patch. By comparing (3.6) and 
(3.7), it can be noted that a stress discontinuity occurs at the interface due to the different 
material properties of the shell and patch. 


(3.7) 


Moments and Forces in the Patches 

By integrating the stresses over the face of a fundamental element, it follows that the 
moment and force resultants for the patches can be expressed as 



with units of moment per unit length and force per unit length, respectively. The explicit 
dependence of the patch’s moment and force resultants on the shell midsurface characteristics 
£x,£o,Kx and Kg can be seen by combining (3.5), (3.7) and (3.8) (we again point out that the 
midsurface strains e x ,£g and curvature changes k x , Kg are unknown and will be determined by 
balancing moments and forces in the combined structure). 


Patch Loadings 

It should be noted that throughout this development, edge effects due to the patches have 
been ignored and thus the expressions in (3.8) apply to patches covering the full circumference 
of the shell and having infinite axial length. The equations can be modified for finite patches 


19 


in the following manner. For a patch with bounding values ij, x 2 , 9\ and 02 as shown in 
Figure 4, the total line moments and forces are 

(M x ) pe = [(Me) pei + (MeU] [^(x) - H 2 (x)] [H,(6) - H 2 {9)) 

(Mt)p' = P,) m + (*/,),„] [II, U) - Hill'll [H,(6) - UAH)] 

( 3 . 9 ) 

(Ay* = ((AT,)*, + (AT,)*,] [H,h:j - H 2 (x)] [//,(«) - UAH)} S K iW)S,;A0) 

(AT,)* = [(AT,)*. + (AT,)*,] [H, (x) - H t (x)] [//,(«) - H t (9)] S,, 2 (x)S,, 2 (0) 

where H is the Heaviside function and H{(x) = H(x — Xi) , i = 1,2, with a similar definition in 
6. The indicator functions Si j2 (x) and Si t2 (0) (see (3.3)) again derive from the property that 
for homogeneous patches having uniform thickness, opposite but equal strains are generated 
about the point (z,0) = ((zi + x 2 )/2, R(9\ + 0 2 )/ 2) in the two coordinate directions. 

The combination of the expressions (3.8) and (3.9) yields the patch moments and forces 
(M x ) pe ,(Mg) pe ,(N x ) pe and (Ng) pe in terms of the middle surface characteristics of the the 
cylindrical shell. Integrating the expressions in (3.8) is somewhat cumbersome however, and 
the procedure can be facilitated by determining the patch moments and forces in terms of the 
resultants of the forced shell. To accomplish this, force and moment balancing is employed. 


Determination of the Patch Moments and Forces 


The application of moment equilibrium about the center of the shell yields the two condi- 
tions 

M x T (M x ) pei T (M x ) P e 2 — 0 

Me + (M e ) pei + (Me) pe2 = 0 (3 ‘ 10) 

where M x and Mq are shell moments. Similarly, force equilibrium in the x and 6 directions 
yields 

N x + {N x )pe 1 + {Nx)p t2 = 0 
No + (N$) pei + {N$)pe 2 == 0 * 

Thus the total patch resultants can be expressed as 


/ h/2 / z \ 

[h/2 

{Mo) P e l + (Me) P e 2 - ~ M e = - / u Q zdz 

J — h/2 

( N x ) pei + (N x ) pej = -N x = - jj a x (l + dz 

[h/2 

( Ng) pei + (N $ ) pe2 = -N e = - / agdz 

J — h/2 


(3.11) 
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which yields 


( A^ar)pei 
{M$) pe i 
pej 

{N$)pe i 


+ (Mx)pe 2 — 

+ (M 0 ) pe2 « 

+ {N x )pe 2 = 

+ ( ) pe 2 ~ 


-Eh 3 

12(1 — // 2 ) 

-Eh 3 

12(1 - i/ 2 ) 

-Eh 


K x + VKg + — J 


, ££.' 
K0 + t'f'-x 


1 — V 2 

-Eh 

1 - V 2 


h 2 k 2 


£x + >'£e + ^ 


h 2 Ke 

£e + U£x + T2R 


( 3 . 12 ) 


The two tangential expressions are approximate in the sense that the terms (1 + z/R) are 
replaced by the truncated geometric series 1 —z/R before integration (this is the same strategy 
which is used when determining the moment and force resultants in the Byrne-Flugge-Lur ye 
general shell theory). The patch resultants in (3.12) are then used in (3.9) to determine the 
total line forces and moments generated by the finite patches. Because these resultants are 
functions of the material properties as well as the midsurface characteristics, they can e 
easily constructed once e x ,£g,K x and k b have been determined. This is again accomplished 
by moment and force balancing. 


Determination of the Midsurface Characteristics 


In terms of the stresses, the moment and force equilibrium equations (3.10) and (3.11) can 
be written as 

CZ (‘ + s) 2,fe + 0 + s) zdz + (' + *) ziz = 0 

rft./2 rh/2+T f-hf 2 

/ (Tgzdz + / {&e)pe l zdz + / {(Tg) pe2 zdz — 0 

J-h/2 Jh/2 J-h/2-T 


(3.13) 


and 


ci °-{ i+! k) dz+ cr^ (‘ + 1) * + o + s) * - ° 


rh/2 rh/2+T f-hf 2 

/ a$dz + / {&e)pe } dz -f / {&o)pe 2 dz — 0 

J-h/2 Jh/2 J-h/2-T 


(3.14) 


rh/2+T 

Jh/2 


The integrals appearing in (3.13) and (3.14) are explicitly evaluated in the [1]. 
After collecting terms, this yields the 4x4 linear system 

[Ashell + Ai + A 2 ] e = Cpeifl + e pe 2 /2 


(3.15) 


in the unknowns e = (e x , e 9 , k x , k 9 ) t . The shell contributions in A sheU , the outer patch 
contributions in A\ and f \ , and the inner patch terms in A 2 and fi are 
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Note that the coefficients o, , i = 1, • ■ • ,4, are of order at most three in h. 
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We point out that A shell , A u A 2 , /i and f 2 depend on material properties of the patches 
and shell (thickness, elastic properties, and Poisson ratios), the radius of curvature of the 
shell, and the voltage being applied to the patches (recall that e pe , = r /31 V\ /T where V 1 is the 
applied voltage into the outer patch with a similar definition for c pe2 ). The above formulation 
isolates the contributions due to the individual patches and is useful if one wants to activate 
only one of the patches. If both patches are present, the above formulation can be simplified 
to yield the linear system 

Ac — Cpe^f] f'2 


where 


i(£l + 2 a 2 ) 
0 


A = 


E 2 + 2 T 


0 

“(El + 2(12) 
E 2 v 4" 2 T Vpe 


E 1 -f 2a 2 
Ei V 4 2a 2 Vpe 



E 2 V + 2 T t/pe E 2 4 2T 


0 


Ei v 4 2a 2 i/ pe 
Ei 4 2,a 2 

0 



Algorithm for Determining the Shell/Patch Interactions 

The steps which are necessary for solving for the shell loads due to the activation of the 
patches can be summarized as follows. 

(1) Set up the 4x4 system At — / and solve for e = (E x ,£e, k x , k$) which contains the 

midsurface strains and changes in curvature. 

(2) Determine the line moments and forces which are generated by the individual patches as 

set up in (3.12). 

( 3 ) The corresponding combined resultants for finite patches are given by (3.9). 

(4) The resultants from (3.9) are directly substituted into (2.13) as the load on the shell if the 

weak shell equations are being used (recall that in this case, q n = 0 and N x = {N x ) pe , 
Ng = (N g )pe, M x = ( M x ) pe , Me = (Me)pe as summarized in (3.1)). For the strong form of 
the equations of motion, the derivative expressions in (3.2) are formed and substituted 
into (2.8) as the external load. 

We point out that the substitution of the patch moments and forces into the strong form 
of the shell equations results in one derivative of the Heaviside and. indicator functions for 
the forces and two derivatives of the Heaviside function for the moments whereas no such 
differentiation is required in the weak form (the derivatives are transferred onto the test 
functions and one simply integrates over the region covered by the patches). This is one 
motivation for using the w r eak form of the shell equations in many applications. 
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4 Patch Contributions to Plate and Beam Equations 

Analysis similar to that used for the thin cylindrical shells can be used to determine the forces 
and moments which are due to the activation of piezoceramic patches which have been bonded 
to a flat plate or beam. 

Plate/Patch Interactions 

The patch interactions with a flat plate can be determined in a manner similar to that 
used in the study of the interactions between a thin cylindrical shell and a pair of piezoceramic 
patches as discussed in the last section. Direct force and moment balancing leads to the 4x4 
system listed under Method 1, and this system can then be solved for the unknowns e x , k x 
and K y . By then substituting these values into resultant expressions similar to those in (3.12), 
one obtains the forces and moments generated by the patch. This procedure can be simplified 
however, by noting that the strains in the x and y directions of a homogeneous flat plate are 
equal when equal free strains are generated by the patch (see [6]). Hence two of the variables 
can be eliminated which leads to the more easily solved 2x2 system given under Method 2. 
It should be noted that the two methods yield the same final force and moment resultants. 

Method 1 

Force and moment balancing similar to that used in the study of the patch/shell interac- 
tions yields the system 


[Ap/ate T Al T A 2 ] e — Cp ei f\ T e p e 2 f 2 

where e = (e r ,£ y , &xi K y) T ■ The matrices and vectors containing contributions due to the plate 
and two patches are 


Aplate 
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e 2 

E 2 i> 


0 E\ E^v 

0 E\v E\ 

E 2 u 0 0 

E 2 0 0 
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&'2 Vpe 
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Tv pe 


® 3 ^pe 
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“ a 3 

-a2Vp e 

«2 

^2^pe 


-a 3 

&2 ^pe 
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T 

T Ope 
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^ 3 Ppe 
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fl 3 Ppe 

-a 3 
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’ «3 " 


' -0-3 ~ 
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-a 3 
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. T _ 
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with 


_ £fe 3 (l ~ i£) 

1 \2E pe (l-v 2 ) 


Ei 


Eh{ 1 - Q 

E pe (l - t/ 2 ) 




Note that these matrices and vectors are identical to those of (3.15) if one takes R — > oo in 
the latter expressions. Once £ x ,£ y ,K x and K y have been determined, the moment and force 
resultants can be found in a manner similar to that used for shells. 


Method 2 

Here we take advantage of the fact that for a homogeneous plate, the strains in the x and y 
directions will be equal when generated by equal free (unconstrained) strains from the patch. 
Hence we take 


t — C x — Cy — 6 + KZ 

as the strain distribution in both the plate and the patch. This yields the stresses 


(J (J x (Jy 


E 

e 

1 - v 


in the plate and 

(^)pe i ~ (^x)pe i (^)pei i ^V e \) 

A Vpe 

( cr )pe 2 = ( cr r)pe 2 = { a y)pe 2 ~~ 1 _ . ( e e P^) 

I f'pe 

in the patch. The force and moment resultants for the patches can be found either by inte- 
grating the stresses over the patches or by using force and moment balancing to express them 
in terms of the resultants for the forced plate. As was done in the shell analysis, we will take 
the latter approach since it yields simpler expressions for the external resultants. Force and 
moment balancing in conjunction with integration of the forced plate resultants then yields 



(ATp)pej T {M x )pe 2 = (A/y) pe] + (^y)pe 2 in/i iA ^ 

( 4 - 1 ) 

, x -Eh 

(A T x )pe l + (A r x)pe 2 = (A^y)pe, + (Ay)pe 2 — j _ • 


As expected, these relations agree with those in (3.12) for the forced shell with the exception 
of the 0(l/R) terms in the latter case which are due to the curvature. 

The total resultants generated by a pair of patches with edges parallel to lines of constant 
x and y can be determined in a manner similar to that used with the shells. For a patch with 
bounding values x x ,x 2 ,y x and 3/2 1 the total resultants are 

(A4V = (M t ) t , = [( + (M.U1 [ff,M - [ H,(y ) - H t (,) ] 

(4.2) 

( N x ) pe = (N y ) pe = [{N x ) pei + (N x ) pe2 } [H x (x) - H 2 (x)] [H x (y) - H 2 (y)} S h2 (x)S 1>2 (y) 
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where again, //,(x) = H(x — x t ) , i = 1 , 2 , Sx, 2 (x) denotes the indicator function described in 
(3.3), and H t (y) and Sx }2 (y) are defined in an analogous manner. 

As before, £ and k must be found in order to determine the resultants in (4.2). This is 
accomplished via moment and force balancing which then yields the system 

[A p iate + Ai + A 2 ] e = e pe ,/i + e p( , 2 / 2 (4.3) 

where e = ( £,k ) T . The component matrices and vectors are 


and 



Ex 

0 



a 2 


a 3 



a Z 
— 0,3 


/. = 


O 3 

T 



The subscripts 1 and 2 again refer to the outer and inner patch contributions, respectively, 
and the constants E \ , E 2 , o 2 and 03 are given by 


E = ~ Vpe) R = Ehjl-U^) 

1 12E pe (l — v) ’ " 2 E pe (\ - u) 


o 2 — 


1 
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-h 2 T + -hT 2 + 
4 2 



a 3 = 


I 

8 



l -T(h + T) . 


Algorithm for Determining the Plate/Patch Interactions 


As in the case of the shell, the steps necessary for determining the plate loads which are 
due to the patches can be summarized in a simple algorithm. 

( 1 ) Set up and solve the 2x2 system Ae = f in (4.3) where e = (e, k) t contains the midsurface 

strain and change in curvature. 

( 2 ) Determine the combined resultants for the finite patches through (4.2) in conjunction 

with (4.1). 


(3) Once determined, the resultants from (4.2) can be substituted directly into the weak form 
of the plate equations (2.15) as the load on the system (with q n = 0 and N x = (N x ) pe , 
N y = (A r y) pe ,M x = (M x ) pe , M y = ( My) pe )• If the strong form of the plate equations is 
being used, the surface loads can be determined via the expressions 


& = -sms,,2W~^ 
* _ 

ni X ^ 


oy 

y ’ 


and these latter values can be substituted into the equilibrium equations (2.14), 
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As in the case of the shells, the use of the strong form results in up to two derivatives of the 
Heaviside function whereas the use of the weak form alleviates this problem by transferring 
the derivatives onto the test functions. 

It should be noted that the voltage choice e pe = e pe , = e pe2 causes pure extension 
(patch pairs excited “in phase”) in the plate while pure bending occurs with the choice 
e pe = — e pe , = e pe2 (“out of phase” excitation). 


Special Case. e p ^ — e P g , — e p ^ 2 

In this case, the constants e and k have the values 

e = 0 


2a 3 

k - — — — e 


fjj -(- 2a 2 


pe 


which leads to the total patch moments 

Th(h + T) 


(A/j,)pe — (A/p)pe 


where 
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This line moment expression is equivalent to the relation 

Pz{ 2 + Pz) 


M X = My 


l+M(3 + p2 + 3p 2 ) 


h 2 7 e pe [H,(x) - H 2 (x)} [H\(y) - H 2 (y )] 


with p z = T/h,h — h/2 which was obtained by Kim and Jones [12] in their development of a 
model for the bending interactions between a flat plate and a piezoelectic. actuator. In their 
work, they consider a patch configuration which excites pure bending in the plate and then 
determine the effective patch moment by first isolating the interface stress of the system. 


Beam/Patch Interactions 


The patch contributions to the dynamics of a thin beam can be determined directly from 
the plate/patch interaction model if one considers only vibrations in the x-direction along 
with the usual transverse vibrations. The system for the beam/patch configuration is then 
identical to that found in (4.3) with the constants E\ and E 2 now given by 


Ex 


h 3 E 

12 


E 2 = h 


E 



Once e and k have been determined, the force and moment resultants for the patch are 
expressed in terms of those of the forced beam and are given by 


(Afr)pei T (A^r)pe2 — El K 
(N x ) pei +(N x ) pe2 = -Eh£ 
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where I — h 3 / 12 is the moment of inertia for the beam. For patches with bounding values x\ 
and x- 2 , the effective moments and forces are 


(Mx) P c = [{M x ) pex + ( M x ) pe2 ] [H(x - x x ) - H(x - x 2 )] 
(N x ) pe = [{N x ) pe] + (iVr)pe 2 ] [H(x -Xi)-H(x- X 2 )]5,, 2 (x) 


which can then be substituted directly into the weak equations (2.18) as loads on the beam 
(with <7„ = 0 and N x = ( N x ) pe , M x = ( M x ) pe ). In order to determine the patch loads for the 
strong form of the beam equations, the corresponding surface moments and forces are found 
via the relationships 


(<7*) pe 


-Si, 2 (x) 


d{N x ) pe 

dx 


( m y)pe 


d(M x ) pe 

dx 


and these latter values are used in (2.17). We again point out that this results in the need to 
differentiate the Heaviside function (once for the force and twice for the moment) whereas this 
problem is avoided in the weak formulation since the derivatives are transferred onto the test 
functions. In fact, the efFect of the Heaviside functions in the latter case is to simply restrict 
the integrals to the region covered by the patches. 


Special Case: Top Patch Activation 


Consider the problem of a beam having only a top activating patch, 
case is 


a 3 

E\ + Cl2 
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03 

e 2 + t 

<*3 


K 

— &pe 
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The system in this 


which implies that the midsurface constants are then given by 


E pe T (E pe T 3 + Eh 3 ) 

£ E 2 e T 4 + 4EE pe T 3 h + 6 EE pe T 2 h 2 + 4 EE pe Th 3 + E 2 h 4 ' Cpei 

6 EE pe Th{T + h ) 

K = EIJ 4 + lEE^THi + QEE^TH 2 + WE^fh 3 + Wh 4 ' ^ ' 

These expressions for e and ac are the same as those found by Gibbs and Fuller [10] when they 
were investigating the moments and forces generated by a single patch which was bonded to 
a thin beam (these expressions are equivalent to their (7) and (8) once the latter have been 
simplified and h/ 2 has been substituted for h). The moment and force resultants (M x ) pe , and 
( Vr) P ei can then be found by substituting e and k into (4.4) (with (M x ) pe2 = (A r r ) pe2 = 0 since 
there is no bottom patch). We note that the resultants in this case are equivalent to those in 
[10] although the forms and signs differ slightly due to a slight difference in the formulation 
of the underlying beam equations. 
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5 Conclusion 


In this work, general models describing the interactions between a pair of piezoceramic patches 
and an underlying elastic structure have been presented. While the presentation is for elastic 
substructures consisting of a thin cylindrical shell, plate and beam, the techniques discussed 
for determining the moments and forces generated by the patches can be directly extended to 
more complex structures and geometries. 

In the case of the shell, the patches are assumed to be curved and the coupling between 
the in-plane strains and the bending, which is due to the curvature, is retained. By using force 
and moment balancing to determine the midsurface strains and changes in curvature of the 
combined structure, expressions for the patch moment and force resultants can be developed. 
In this manner, the loading due to an actuating pair of patches can be expressed in terms 
of the material properties of the shell and patches (thickness, elastic properties and Poisson 
ratios), the radius of curvature of the shell, and the voltage being applied to each of the 
patches. This provides a shell/patch interaction model which retains the curvature effects 
as well as admits potentially different voltages into the two patches. We point out that the 
general techniques used for determining this cylindrical shell/patch model can also be used 
to determine the interactions between pairs of piezoceramic patches and more general shells 
(for example, in the case of a spherical shell, one would retain the curvature effects in both 
coordinate directions). 

The techniques for determining the patch interactions with a cylindrical shell were then 
used to develop general interaction models for patches which are bonded to thin flat plates 
and beams. As in the shell case, the models are sufficiently general to allow for potentially 
differing patch voltages which implies that they can be used for controlling system dynamics 
when both flexural and extensional vibrations are present. To compare with existing analyses, 
the plate/patch interaction model is shown to be equivalent to that of [12] in the special case 
when pure bending motion is excited. Also, the beam/patch model is equivalent to that of 
[10] when there is one actuating patch. Hence the beam and plate interaction models are 
consistent with existing theories in the special cases previously examined while also allowing 
for more general structure/patch interactions which can arise in more complex applications 
(for example, coupled systems). 

For each of the shell, plate and beam interaction models, the contributions of the patches 
are carefully described in both the strong and weak forms of the time-dependent structural 
equations of motion. This provides models which can be used in a variety of applications 
including numerical simulations, parameter identification, and control schemes. In each of 
these applications, the models are sufficiently general to provide for a variety of approximation 
techniques including modal, spectral, spline and finite element schemes. Finally, the patch 
loads determined by these interaction models can be applied to higher order structural models 
in exactly the same manner, and analogous models can be used for multiple patch pairs and 
more complex geometries. 
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